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ABSTRACT 

An experimental investigation was made to determine the aerodynamic 
response of a wind anemometer to unsteady air loads. The anemometer, a 
perforated drag sphere mounted upon a cylindrical pedestal, sensed 
dynamic pressure by a strain gage balance system. 

The anemometer response was calibrated for static and dynamic 
mechanical inputs and for steady aerodynamic loads prior to the un- 
steady investigation. The unsteady environment was obtained in an 
oscillating flow wind tunnel 1 iioiaes a rotating shutter valve to 
superimpose a harmonically varying velocity upon a mean free stream. 

For steady flow conditions the anemometer drag coefficient remained 
essentially constant, apparently independent of Remiolds number. The 
anemometer responded well to the unsteady air loads at frequencies up 
to 20 cycles per second and mean dynamic pressures to 8.85 pounds per 
square foot. Unsteady aerodynamic response was qualitatively substan- 


tiated by a quasi-steady approximation to the measured unsteady drag. 
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ERRATA 
14 eee eesomould read contrast vice constrast. 
17, line 26: should read shown vice fhown. 


106, Lett Nand sketch: should read steel 
rod vice steel rid. 


JO09, bottom line: adda period at the end 
of the sentence. 


110, line 3 of the 3rd paragraph: the symbol 
w should be inserted in order to read "...the 
prOodtg@iwwes is...". 


1137 eet ine Under Equation E.04: 4should 
read E.04 vice E,04. 
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1l.. MSmieroduction 

Structures exposed to ground winds can experience, depending on 
the geometry and/or motion involved, both steady and unsteady bending, 
torsional, and distortional loads. The effects of wind loads upon the 
elastic response of large structures have been observed as the oscilla- 
tions of tall slender smoke stacks and antennas, "galloping" of trans- 
mission lines, warping motion of suspension bridges, and more recently 
as the vibrations of space launch vehicles while on the launch pad. 
These effects have been the subjects of such investigators as Den 
Hartog [ 3 ] » Scruton and Walshe [10] » scruton (11) and Buell, et at [1]. 

Many of the studies conducted on ground wind loads have used wind 
tunnel tests in which the velocity was steady in magnitude and direction 
and whose profile was essentially uniform (tunnel boundary layer was 
small compared to model dimensions). In reality the ground wind is 
rarely constant in either magnitude or direction and the boundary layer 
of a ground wind can be several hundred feet thick. A structure's re- 
sponse to a steady uniform flow can be quite different from its response 
to a flow which varies in space and time. Although correct reproduction 
of aerodynamic flow properties as a prerequisite for meaningful dynamic 
model tests has been recognized for some time, possibly the first tests 
in which the dynamics of the structure and the turbulent boundary layer 
flow were both reproduced occurred in 1964, directed by A. G. Davenport 
fe2 | in cooperation with J. E. Cermak using the Army meteorological wind 
tunnel at the Colorado State University. 

Any meaningful prediction of wind loads must be based on a know- 
ledge of the wind itself, which in turn implies accurate measurements 
of the wind. These measurements are becoming increasingly more important 


as the studies of wind loading on large elastic structures become more 
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detailed. Measurement of fluctuating, gusty winds is difficult with 
instruments that are frequency-response limited, and hence can only 
indicate averages of the various components. Recently the National 
Aeronautics and Space Administration, Langley Research Center (NASA 
Langley), Hampton, Virginia, has developed a fast response anemometer 
capable of following wind fluctuations in the range of zero to five 
cycles per second. According to Reed and Lynch [9] this development 
resulted from the gratifying performance of a ping-pong ball, attached 
to a single component force pick-up, which was used to measure time 
fluctuations of low level dynamic pressures in the early 1950's. Ad- 
ditionally, wind tunnel testing at NASA/Langley demonstrated that the 
Steady drag force vector for a perforated hollow sphere remained 
aligned with the steady wind velocity vector within a tolerance of +e 
degrees. In constrast a smooth sphere exhibited an abrupt, random 
shift of the drag vector to either of two poeterons . symmetrically 
located about 15 degrees from the wind direction. The latter result 
is not surprising since the separation point for the flow is not geo- 
metrically stabilized as is the case for a perforated sphere. The choice 
of a perforated sphere for use in a drag anemometer is a logical one 
based on these observations. 

The drag sphere anemometer, in addition to fast response, has 
several other important characteristics: 

(1) measurable drag forces are generated in the low dynamic 

pressure regions. 
(2) the coefficient of drag is constant over a wide range of 
velocities, hence dynamic pressure can be readily interpreted 


from the drag. 
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(3) the natural frequency is high compared to the frequencies of 

interest. 

(4) the drag force vector remains aligned with the wind direction. 

Current use of the drag sphere anemometer has been indicated by 
Duncan and Foughner [ 4 ] as part of a continuing ground wind loads in- 
vestigation at Wallops Island, Virginia. 

An oscillating flow may be defined to consist of a harmonically 
varying velocity superimposed upon a constant free stream. Low speed 
oscillating flow investigations have, in general, been limited. Miller 
fie? has investigated boundary layer transition over a flat plate. 

Murphy [ 8 | recently investigated the pressure distribution about a 
circular cylinder in an unsteady flow. No investigations concerning 
perforated spheres in an oscillating flow were found in the literature. 

This investigation was undertaken to determine the aerodynamic 
behavior of the perforated-sphere in an oscillating flow. The frequency 
range of interest was from zero to twenty cycles per second, with dynamic 
pressures to 8.85 pounds per square foot, which approximately corresponds 


to air velocities of 88 feet per second, or 60 miles per hour. 


ae, Experimental Equipment and Procedures 


The perforated sphere anemometer used in these investigations was 
supplied to the Naval Postgraduate School by the Aeroelastic Branch of 
NASA (Langley) in support of the aerodynamic frequency response program. 
The anemometer was floor mounted when installed in either the 3.5 x 5.0 
foot low-speed wind tunnel or the 2.0 x 2.0 foot nonsteady wind tunnel, 
both of which are located at the Naval Postgraduate School, Monterey, 


California. 
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The seven-inch drag sphere was made of plastic (weight and corrosion 
considerations) and had forty 0.75-inch diameter perforations uniformly 
distributed over the surface (see Figure 1). It was mounted on a two- 
component strain gage balance system oriented such that the balances 
sensed the two horizontal components of sphere drag force. Static in- 
ertial compensation was provided by an internal counterweight arrange- 
ment which could be positioned such that the center of gravity of the 
sphere system corresponded to the effective center of the strain gage 
balances. No adjustments of the balance weight were made Een owminn ¢ 
receipt of the anemometer from NASA/Langley. 

Anemometer Calibrations 

The anemometer strain gage balances, conventional 350 ohm full- 
bridge circuits, were statically calibrated through force ranges that 
corresponded to gage strain levels of approximately + 1000 micro-inches 
per inch. With a bridge excitation level of 10.0 volts (D.C.), the 
sensitivities were found to be approximately 4.42 and 4.15 millivolts 
per pound respectively, with the results being linear in the range in- 
vestigated. See Figure 2 for actual calibration curves. 

Although both balances were calibrated in order to allow simul- 
taneous evaluation of lateral (side) and drag forces, only one balance 
was operative following the calibrations, since the second balance be- 
came unusable due to an apparent unbonding of one of the strain gage 
elements. Consequently the scope of the program became restricted to 
a frequency response evallation of sphere drag force. 

Astrodata differential D.C. amplifiers, model A885, were used in 
these calibrations, with a gain of 100 to increase the output level to 


a value suitable for magnetic tape recording, and to allow biasing of 
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the steady drag component so that only the unsteady portion of the signal 
would be recorded. Balance sensitivities stated in a previous para- 
graph are equivalent to amplifier input values. 

Since the anemometer was to be exposed to velocity perturbations 
while attached to the floor of the nonsteady tunnel, which definitely is 
a vibratory environment, it was necessary to consider the inertial 
response traits of the anemometer for possible correctional terms 
during data reduction. The need for such considerations can be visual- 
ized by realizing that the strain gage balance is an elastic system 
that relies upon:relative motion between the sphere and the pedestal 
base in order to produce measurable strains in the balance system. 

It would normally be necessary to consider both rotational and 
translational inertia terms, however, during preliminary tests the non- 
steady wind tunnel (described later in this section) was found to have 
floor motion predominantly in a direction tangential to the floor with 
maximum accelerations of 2.0 "g's" at 10.0 and 17.5 cps. This accelera- 
tion direction is dominant since the tunnel test section is relatively 
stiff, and the tunnel is mounted on a framework of steel pipes which 
act like a "bent", without significant bending moment constraint at 
the joints. Therefore a calibration was made only in the translational 
vibration mode. 

The inertial frequency response calibration was performed over a 
frequency range of 5 to 100 cps. A Calidyne Shaker, model 219, operated 
by a Ling control console, model CP 3/4 was used to accelerate a steel 
plate upon which the anemometer and reference accelerometer were mounted, 
as fhown in Figure 1. Horizontal accelerations were measured with a 
Statham accelerometer, model SCO 210, which was used as a reference 


sensor unit in both the mechanics laboratory and the nonsteady wind 
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tunnel. The primary range of frequencies was below 100 cps, and this 
was readily measured by the Statham strain gage device since it had 
a uSable frequency range from O(D.C.) to 150 cps, and an acceleration 
range of + 5 "g's" (see Appendix A for accelerometer calibration). The 
shaker was connected to the plate by a 1/16 inch diameter steel wire 
flexure coupling, which had an estimated critical buckling load of 
5400 pounds. A wire flexure coupling was used to avoid imposing trans- 
verse loads onto the shaker table. The steel plate rested upon a steel 
surface table which was adjusted to be level. The surface table and 
the plate had polished surfaces, and an oil film was applied between 
the moving surfaces in order to reduce the level of frictional forces. 
The coupling was positioned such that the axial force transmitted 
through the wire coupling acted approximately through the center of 
mass of the combined dynamic assembly (plate, anemometer, and accelero- 
meter). 

If the anemometer were represented as a system with a single degree 
of freedom: ai 


K 


YY) BASE 


a | 
x(t) —» yt) 


then, the equation of motion for the freely vibrating mass is: 
mx + C(k-y) +k(x-¥)= ©. (2.01) 
Loewen , and rewrite 


me+ C2+kz = -™MY 
then for sinusoidal motion of the base, ar Yo sinvwot 


z me 2 Fy =e Wiz = Yo Ly® Sin wt (2.02) 


Ls 


where 


202 , we 

The solution of (2.02) is: 

z(t)= cata Sin(wt- d)] . 
x 7% La) \% 
wh I -(SY]+ 49S) 


D — arctan a a 


- Gy’) 


where 7- (t) is the displacement of the anemometer force pickup rela- 


— 


tive to the base, and Yow is the magnitude of the harmonic forcing 
function. 

Accelerations of various magnitudes and constant frequency were 
applied to the base plate at four intervals up to and down from the 
maximum of 2.0 "g's''. The anemometer response 4 (t), accelerometer 
output, Yo, and phase difference between these sensors were recorded 
at each level of acceleration. The frequency was increased and the 
procedure repeated. Figure 3 is a summary of the anemometer inertial 
response calibration. 

The anemometer natural frequency was found to be approximately 
60 cps as indicated by maximum response from oscilloscope traces, and 
by observation of large excursions of the perforated sphere. For 
frequencies greater than 60 cps, the response wave forms became very 
distorted. Considerable difficulty was encountered in the determination 
of a 90 degree phase shift between input and output. The difficulty 


was attributed to the anemometer not behaving as a system with a single 


degree of freedom. 
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Anemometer damping was evaluated by observing the output response 
to a step input. No low frequency oscillation ensued for any of the 
step inputs, and it was concluded that the system fundamental mode was 
at least critically damped. Examination of the response, however, with 
the single sweep capability of the Tektronix (type 551) oscilloscope, 
revealed a high frequency (greater than 100 cps) signal superimposed 
upon the basic wave form. The excitation of this high frequency signal 
was controllable to a degree and was apparently caused by a ringing 
mode of the perforated shell. i 

Nonsteady Tunnel Investigation 

The investigation of anemometer unsteady aerodynamic response was 
conducted in the low speed, oscillating flow tunnel located in the 
Aeronautical Laboratories of the Naval Postgraduate School. The open 
circuit tunnel has a 24 inch square test section, 223.5 inches long. 
Air is driven by two series-connected Joy Axivane fans connected 
directly to 100 hp, 1750 rpm electric motors. The tunnel has a bell- 
mouthed entrance with a nozzle area contraction ratio of 16:1. There 
are three fine mesh screens in the settling section area which reduce 
steady flow turbulence levels in the test section below 0.1 percent 
of mean velocity. 

Construction of the test section floor and ceiling is from 2.0 inch 
thick aluminum, the sides from 2.0 inch lucite. Heavy construction was 
necessary to reduce structural deflections caused by large variations 
in static pressure. 

Oscillating flow was created by a rotating shutter valve assembly 
located at the trailing edge of the test section and driven through 
belts by a variable speed electric motor. The shutter blades were 


removable to permit installation of various widths that could enclose 


20 


from 50 to 100 percent of the test section area, and which provided 
unsteady velocity amplitudes from 1.3 to 84.0 percent of free stream 
velocity. Various pulley ratios of the belt drive were available 
effecting an operational range of frequency from 2 cps to 933 cps. 

Murphy [ 8] performed steady flow velocity surveys across the test 
section and found that the variation of V relative to van was less 
than + 1 percent for distances greater than three inches from the 
tunnel walls. 

For the actual nonsteady tunnel tests it was desired to subject 
the aneometer to five values of mean dynamic pressure with each of the 
three available shutter vane sizes. The investigation in general in- 
cluded fifteen points of data for each run at various frequencies up 
te 20) cps. 

In order to determine the effect of blockage produced by the ro- 
tating shutter vanes, and to insure that the hot wire system would per- 
form satisfactorily, several trial runs were performed with the anemon- 
eter absent from the test section. Dynamic pressures of 7.8 and 15.6 
pounds per square foot (psf) were set with the vanes horizontal. The 
hot wire DC output was set to read full scale (1.0 volts) for zero 
frequency and maximum velocity. The tunnel was operated at various 
frequencies and the drop of the DC level of the hot wire was noted. 
Figure 4 is typical of results for an initial mean dynamic pressure of 
7.8 psf. Post-run checks were performed to ensure the absence of hot 
wire drift. 

With the information obtained from an evaluation of the blockage 
effect, a zero frequency dynamic pressure was selected. The hot wire 
DC attenuation was adjusted to full scale, and the zero frequency RMS 


values were recorded for the anemometer, hot wire, and accelerometer. 
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The speed of the rotating shutter valve was initially adjusted to the 
minimum, 3.6 cps, for the combination of pulley ratios used. After 
amplified anemometer and hot wire DC values were biased to zero and 
the oscillation frequency was stabilized, tape recording of signals 
for the particular data point was commenced. RMS values of the ane- 
mometer, hot wire, frequency, and mean dynamic pressure were logged 
manually. This procedure was repeated at seven increments of frequency 
up to and down from the maximum of 20 cps. Post-run zero frequency and 
zero dynamic pressure were noted. This procedure of data acquisition 
was repeated at five levels of dynamic pressure for each of the shutter 
vane widths. 

Nonsteady Tunnel Instrumentation 

The Ampex magnetic tape system used was an electron tube type 
model FR-100 with 14-channel recording capability. Both frequency 
modulated and direct record amplifiers were available. A tape speed 
of 7.5 inches per second was used throughout the data acquisition, 
which provided a flat frequency response on the FM channels from 0 to 
1250 cps. Three frequency modulated channels were used to record ane- 
mometer, hot wire, and accelerometer signals, whereas voice information, 
to identify blocks of taped data, was recorded through a direct record 
amplifier. 

A Tektronix type 555 oscilloscope with a pair of 4-channel pre- 
amplifiers was used to simultaneously monitor original and recorded 
data. 

An integrating digital voltmeter (Vidar 510) was used to alter- 
nately monitor the DC output of the anemometer bridge circuit and the 


hot wire anemometer. The DC values of the amplified signals from both 
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instruments were maintained at values as near to zero as possible by a 
battery-supplied bias voltage. The DC biasing permitted the tape re- 
cording of optimum level AC signals. Figure 5 presents a block diagram 
of the entire data acquisition system shown by photograph in Figure 6. 
Details of the anemometer installation are shown in Figure 7. Shielded 
co-axial cables with BNC fittings were used to reduce the influence of 
extraneous signals, for example 60-cycle ground loops. A centrally 
located switch panel afforded an orderly flow of data and convenient 
monitoring of information. 

The hot wire probe was constructed at the Postgraduate School by 
soldering an 0.0015 inch diameter tungsten wire on a pair of jewelers 
broaches. The jewelers broaches were i ctpegata in epoxy plastic, which 
was fitted to a 1/4 inch hollow steel rod. The performance of the 
Security Associates bridge circuit/analog computer system was excellent. 
Repeated results assured a linear DC variation from zero to full scale 
and a high degree of system stability was noted. The sensitivity of 
the system was such that outside wind fluctuations during steady opera- 


tion of the wind tunnel were detectable. 
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3. Experimental Results and Discussion: 


Run Program 
The investigation of the aerodynamic response of the wind anemometer 
performed in the oscillating flow tunnel was separated into three phases: 
(a,) tunnel familiarization and calibration of flow, (b.) evaluation 
of the steady flow drag, and (c.) the determination of unsteady drag 


in oscillating flow. The phases may be further subdivided as follows: 


Phase I Calibration 

Run q(psf) f(cps) Objectives ‘ 

Cl fro” = Oe ae Tunnel familiarization, check blockage 
effects on small vanes, and evaluate hot 
wire performance 

C2 - Medium vane blockage 

C3 tt tt Large tt it 

C4 Teo , Smal” 

C5 if tt Medium tt tt 

C6 tt tt Large tt tt 

Phase IL Steady Drag 

Run q(psf) f (cps) Objectives 

C7 Poel os 7 0 Determine anemometer output vs q 

cs 1.0-15.6 Ut tt it tt UL 1 

C9 Zee, 0) VOv=" 19.0 Determine velocity perturbation effects 
upon mean value of drag 

Phase III Unsteady Drag 
Run q(psf) f (cps) Ob jective 
1-4, 1.0-7.6 3.6 - 19.3 Effects of small vanes (small velocity 
10 perturbation) upon unsteady drag 

5-9 U0.6-1.7 3.6-- 21.0 Effects of large vanes upon unsteady drag 
“Teer 0.7-2.9 i i tmeatum i a 4 
16-18 4.2-7.8 3.6 - 19.1 a "small " - 7 . 
ZO eee x '" medium " e u o 
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Data Reduction 

The initial step in the reduction of experimental data was to 
evaluate the operating conditions for which the anemometer output was 
Significantly affected by tunnel floor accelerations. In consideration 
of the inherent fluctuations of anemometer RMS output, "significantly" 
was interpreted to mean those data points which influenced the total 
anemometer signal by more than five percent. Frequencies which caused 
excessive tunnel vibration were avoided, but using the criteria above 
it was necessary to correct about 25 percent of the data points. 

A method to correct the anemometer response, both in magnitude 
and phase, is presented in Appendix B. An important assumption in the 
analysis is that the dominant floor accelerations are of a sinusoidal 
nature and of the same frequency as the oscillating air flow. The 
assumption was based upon oscilloscope photographs that were made of 
the response signals reproduced from the tape recorded data. Two Kron- 
Hite band pass filters were employed in the playback channels to elim- 
inate the high frequency components of the actual signals. Tables l, 
2 and 3 contain the corrected data for small, medium and large shutter 
vanes respectively. 

Appendix C is a development of hot wire response to a base or 
floor acceleration based upon elementary beam theory. For all values 
of frequency and acceleration in the range investigated this effect 
can be neglected. The effect of tunnel floor velocity on hot wire 
response, considered in Appendix D, is less than two percent of the 
unsteady air velocity and therefore negligible. The hot wire output 
can then be considered as the true behavior of the oscillating air 


flow. 
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The phase differences between (a.) anemometer and hot wire signals 
and (b.) anemometer and accelerometer signals were determined by 
averaging the output of an AD-YU phase meter for an interval of taped 
data (30-45 secs). The difference was qualitatively substantiated by 
photographic measurements. Phase angles of the anemometer relative to 
the hot wire are listed in Tables 4 and 5 for the small and medium 
shutter vanes. Phase differences for the large vanes are omitted due 
to the large degree of variation resulting from the complex waveform. 

Drag-Steady Values : 

Figures 8 and 9 are plots of raw anemometer drag as a function of 
dynamic pressure as obtained in the steady West Coast Research (W.C.R.) 
tunnel and in the oscillating flow tunnel for zero frequency. The 
steady coefficient of drag was evaluated for each set of data using 
the anemometer static sensitivity (Figure 2) previously determined. 

The drag coefficient obtained from the W.C.R. tunnel was 0.69, while 

the non-steady tunnel yielded 0.71, which includes an estimated cor- 
rection for solid blocking. For both evaluations, the anemometer out- 
put was a linear function of dynamic pressure in the range investigated, 
which indicated that the drag coefficient for the perforated drag sphere 
was apparently independent of Reynolds Number. 

Reed [9] has shown for the one-dimensional case of a non-steady 
velocity v(t) superimposed upon a mean free stream V that the time 


average of the instantaneous velocity squared is 


Vis Vi(I+ 5) 
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Taking the square root of both sides, expanding the right hand 

side in a binomial series, and retaining first order terms yields: 
VF evore Ss] 

In other words, the mean velocity when determined by a measurement of 

dynamic pressure (or v2) is biased by the turbulence in the flow. 

The unsteady tunnel was used to determine the effect of an oscil- 
lating flow on the mean drag of the anemometer. The first step was to 
measure the drag in a steady flow over a dynamic pressure range of 
1.04 to 15.6 pounds per square foot. This test corroborated earlier 
results obtained in the W.C.R. tunnel. Next the mean anemometer drag 
was recorded in an oscillating flow (both small and medium sized vanes 
were used) over a mean dynamic pressure range of 2.08 to 8.84 pounds 
per square foot. The mean drag was recorded at about 16 discrete 
frequencies between 3.6 and 19.1 cycles per second at each of 9 dynamic 
pressure levels. It was found that the mean anemometer drag at a mean 
dynamic pressure was independent of flow oscillation frequency and 
equal to the steady flow value. The results of the steady and the un- 
steady flow tests are shown in Figure 9. 

Unsteady Drag and Phasing 

The variation of unsteady drag as a function of frequency and 
dynamic pressure for the case of the small shutter vanes is presented 
in Figure 10. The individual data points for the curves were obtained 
by positive incremental frequency steps up to the maximum obtainable, 
and then by negative steps to the minimum frequency. A high degree of 
repeatability was obtained for the small vanes. Unsteady drag, for 
the lowest dynamic pressure (1.0 psf) in the series, remained essentially 


constant with frequency, while each of the other curves for higher 


ai 


dynamic pressures generally decreased with frequency and apparently 
approached a limiting value. It is interesting to note that the varia- 
tion of unsteady drag with dynamic pressure at a constant frequency 
of 19 cps is approximately one-fifth the value at a frequency of 
3.6 cps. 

All of the curves exhibited a peak in the vicinity of 10 cps, 
the relative magnitude of which increased with dynamic pressure. Also, 
at frequencies near 10 cps, the ratio of V/V was a local maximum for 
the particular run, and is believed to cause the rise in unsteady drag. 
The reason for an increase in this velocity perturbation ratio at this 
frequency is not known but was attributed to the flow characteristics 
of the tunnel. 

The unsteady drag variation for the medium and large shutter 
vanes is shown in Figures 11 and 12. For the case of medium vanes 
an additional peak appears at approximately 17.5 cps, which also can 
be attributed to a local maximum of ie Although the medium-vane 
curves display a lower degree of consistency than do the small-vane 
curves, the trend of decreasing magnitudes with frequency is maintained, 
as is the trend toward a limiting value for each dynamic pressure with 
increasing frequency. The ee shutter-vane results display an even 
lower degree of consistency, but still display the drag peaks at 
approximately 10 and 17 cps. Unsteady drag dependence on dynamic 
pressure, as indicated by the smaller vane data, is substantiated for 
frequencies less than 11 cps; for frequencies greater than 11 cps the 
scatter of data precludes immediate interpretation. 

The anemometer phase lag to an aerodynamic forcing function is 
affected both by dynamic pressure and frequency as shown in Figure 13, 


which is a summary of results for the small velocity perturbations 
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(V/V = 0.02). Phase lag generally increases with frequency and 
apparently has no dependence upon dynamic pressure. The increase in 
phase lag with frequency up to about 16 cps was also observed during 
the inertial investigation described in Section 2 (cf,. Figure 3). 

With the exception of the lowest dynamic pressure (0.73 psf) in 
the medium vane runs (where 0/V = 0.12), anemometer phase lag re- 
mained near zero up to a frequency of 12 cps, where a rapid increase 
was observed toward a peak value occurring at approximately 15 cps. 
See Figure 14. 

Quasi-Steady Assumption 

A theoretical treatment of bluff body oscillating flows is very 
scarce in published literature. Isaacs[ 6] ,» in an investigation of 
the non-steady flow over a helicopter blade at small angle of attack, 
has approached the problem for a potential flow case. However, the 
analysis included two limiting assumptions: the flow was two-dimen- 
sional over a flat plate airfoil, and the vortices shed remain in the 
plane of the airfoil. The method is not directly applicable to the 
three-dimensional sphere problem because of obvious differences be- 
tween idealized potential flow and bluff body separated flow, but 
Isaacs does suggest that the lift on the airfoil can be evaluated 
from an averaging of the instantaneous lift due to the variable velocity. 

Fung [ 5 ] further defines this quasi-steady principle, which assumes 
that the characteristics of an airfoil in an unsteady flow at any in- 
stant of time are equal to the characteristics of the same airfoil in 
a steady flow with flow parameters equal to the instantaneous values. 

Appendix E is a development of the quasi-steady principle for the 


three-dimensional sphere. 


yas 


Drag Ratio 

Anemometer drag ratio (DR,) is defined as the quotient of un- 
steady drag divided by drag produced by the steady component of velocity. 
The quasi-steady drag ratio (DR) is approximately equal to 7 and is 
developed in Appendix E—. These ratios are plotted versus Strouhal 
number (A+) to indicate the variation with frequency and mean velocity, 
and to evaluate the degree of validity of the quasi-steady assumption. 
The scales for the coordinate axes are not consistent due to the large 
variation of maximum values between runs. - 

The quasi-steady assumption agrees best with the anemometer 
response for the runs with the small shutter vanes installed. Two 
factors are significant for this result: (1) velocity perturbations 
produced by the small vanes are better approximations to harmonic 
functions, and (2) the ratio of v/v is small, which is inherent to 
the approximation of quasi-steady theory. Photographs of response 
signals from oscilloscope traces are presented in Figures 15 and 16 to 
illustrate the effects of shutter vane size upon the deviation from 
harmonic perturbations. 

The variation of drag ratio with Strouhal number and dynamic 
pressure for five runs with the small shutter vanes installed is con- 
tained in a sequence of figures beginning with Figure 18. The ane- 
mometer response curves have the same form as the curves for the quasi- 
Steady approach (based upon the hot wire output) or, in other words, 
the anemometer is sensing the velocity perturbations throughout the 


frequency range. The level of response increases, as expected, with 


increasing dynamic pressure. However, for a given dynamic pressure, 
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both the average value of the drag ratio and the quotient, DRy/DRo, 
decrease with increasing Strouhal number. Figures 20, 21 and 22 also 
contain results obtained from separate test dates to demonstrate the 
repeatability of data and thereby provide a confirmation of results. 

Drag ratio variation obtained with the medium vanes installed is 
presented in a second sequence of figures, starting with Figure 23. 
The correlation of anemometer and quasi-steady results is still main- 
tained, although the appearance of a second peak occurs at a Strouhal 
number of 0.22 (the equivalent frequency is approximately 17.5 cps). 
For the medium vanes, the average value of the drag ratios are on the 
order of seven times the small vane values, although the actual vane 
size has been increased only 50 percent. 

Large vane effects upon drag ratio exhibit more deviation and 
less repeatability than the vane sizes previously considered. The 
large vanes create more distortion. of the velocity perturbation, and 
anemometer response becomes more random in nature. Large fluctuations 
(50% of the freestream value) produced anemometer response that was 
generally about half of the quasi-steady value, as shown in the series 
of figures starting with Figure 28. 

Due to the non-linear variation of the tunnel nonsteady velocity 
for a given mean dynamic pressure and shutter size, as a function of 
frequency, it is not clear what the effects are, if any, of aerodynamic 
damping upon anemometer response. However, a tendency toward reduced 
response (as evidenced by the reduced drag ratio in relation to the hot 
wire value) has been observed for a combination of the greater dynamic 


pressures, large fluctuations, and higher frequencies. 
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Throughout the nonsteady investigation, the tunnel radiated 
various levels of noise dependent upon frequency, shutter size, and 
dynamic pressure. If the inlet test section is compared to an open 
organ pipe, then the fundamental frequency may be calculated equal to 
27 cycles per second. The peaks of the drag ratio curves may then re- 
flect a significant degree of sub-harmonic response. 

Harmonic Considerations 

The output signals obtained from the anemometer and hot wire were 
obviously not true sinusoids. Observations of these signals indicated 
a possible significant higher harmonic component, as shown in Figure 33. 
The anemometer response has been shown to follow the gross aerodynamic 
input reasonably well but a more quantitative definition of response 
was desired. Therefore, the gross signals were examined to determine 
the amplitude and phase relationships among the fundamental, second, 
and third harmonic components and between the anemometer and hot wire 
Signals. A numerical approximation of a Fourier analysis was used, as 
indicated in Appendix F. Since the required computations were performed 
with a desk calculator, only a few selected data points were analyzed. 
Results of this examination are shown in Table 6. As a check on the 
consistency of the approximations, and the results, two different data 
samples from runs 1-2 and 2-4 were analyzed. The data samples from 
run 1-2 were taken from different cycles on the same oscilloscope 
photograph, while the samples from run 2-4 were taken from separate 
photographs. Quantities with good correlation, and discrepancies, are 
obvious in Table 6. The eocticie of a significant second harmonic 
component is definitely shown, having an average amplitude, over the 
points tested, of 23 percent and 17 percent of the fundamental component 


amplitudes for the hot wire and anemometer respectively. Corresponding 
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values for the third harmonic component are 4 and 6 percent. The 
fundamental component of the hot wire led that of the anemometer by 
an average of 9 degrees. The phase angles between hot wire and ane- 
mometer computed for the second and third harmonic components were 
random in sign and magnitude. 

Inspection of several oscilloscope photographs (cf. Figure 34) 
revealed a higher frequency superimposed on the basic signal. This 
higher frequency, about 145 cps, shows up on both hot wire and ane- 
mometer signals, but it is much more prominent for the latter. The 
sensitivity of the anemometer to this input frequency is believed due 
to the oscillations of the sphere itself, as opposed to the sphere- 
pedestal combination. This "ringing'' of the sphere was also noted 
during part of the inertial response calibration to determine ane- 


mometer damping ratio. 


33 


4. Conclusions 
The natural frequency of the anemometer was 60 cps with approxi- 
mately a critical amount of damping. 

For steady flow conditions the anemometer has a constant drag 
coefficient, thus apparently independent of Reynolds Number within 
the dynamic pressure range of 0.78 to 15.6 psf. 

The average drag measured by the anemometer in an unsteady flow 
for a given mean dynamic pressure was independent of frequency and 
equal to the drag measured in a steady flow at the same dynamic 
pressure. 

Unsteady drag, as determined from a quasi-steady approach, pro- 
vided a good qualitative comparison to the unsteady drag measured 
by the anemometer. Quantitative results of the quasi-steady assumption 
provided only an order of magnitude indication of anemometer response. 

The unsteady component of anemometer response increased with mean 
dynamic pressure and with unsteady velocity. For the runs performed 
with the small shutter vanes, the unsteady anemometer response rela- 
tive to the unsteady velocity decreased with increasing frequency, 
and apparently approached a limiting value. 

The investigation was complicated by an amplification of the un- 
steady air flow at several frequencies in the range of interest, and 
by a significant second harmonic component, about 25 percent of the 
fundamental amplitude. 

It is recommended that future investigation be directed toward a 
thorough definition of the oscillating flow behavior of the nonsteady 


tunnel throughout its operating range. 
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Run* 


1-1 


I= 


1-4 


less 


1-6 


ey/ 


1-8 


q 
(CPS) (psf) 
3.6 4.21 
5.0 dy 2a 
7 4.21 
10.0 Gaal 
13.0 4. Dil 
lige 8 4.21 
19.2 d, Dal 
17.8 4.21 
15), 2 4.16 
11.6 4.16 
9.2 4.16 
6.4 4.16 
4.2 Ly oG 
30 4.16 


*The run number 


TABLE 1 


REDUCED DATA: SMALL VANES 


is followed by a coordination number which 


(VRMS) 
x 102 


0. 


0. 


OF 


QO. 


CA 


031 


026 


.018 


822 


.013 


Olt 


.010 


O11 


012 


.016 


shu, 


020 


029 


031 


to identify data as recorded on magnetic tape. 
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<\l<? 


068 
054 
042 
056 
039 
.036 
044 
039 
036 
052 
048 
046 
062 


.070 


is used 


TABLE 1 (Cont'd) 


v 
Run ; : " “a 
(CPS) (psf) Ueis ) 
scauliO 
2-1 omG 2.34 0.0130 0.044 
2-2 6.0 oe 0.0100 0.032 
2-3 8.8 2.34 0.0084 0.031 
2-4 Mlege2 2 3G 0.0090 0.034 
25 323 Ze ot 0.0068 0.026 
2-6 LO 2.34 0.0068 0.025 
2-7 16.9 ona 0.0064 0.075 
2-8 1 9p) Zeeoe 0.0064 0.4028 
2-9 onl 2.34 0.0060 0.027 
2-10 159 2.34 0.0065 0.025 
aod. 14.1 2.34 0.0066 0.025 
2-12 10.1 ae 0.0099 ~ 091089 
2-13 70 a o4 0.0090 0.028 
2-14 5.0 DOG 0.0110 0.034 
2-15 4.2 2.34 0.0120 0.038 
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TABLE 1 (Cont'd) 


<l<?7 


Run 
(CPS) (psf) (VRMS ) 


x 102 


coe anc 7 160 0.077 0.087 
oo 5.8 7.60 0.050 0.064 
323 8.6 7.55 0.037 0.056 
3-4 wih 7.55 0.047 0.070 
3-5 1806 7.55 0.026 0.054 
$26 17ea 7.55 0.017 0.058 
327 19.2 7.60 0.018 0.078 
3-8 15.7 # 60 0.019 0.054 
3-9 12.4 7.55 0.033 0.058 
3-10 9.7 7.55 0.050 0.076 
3-8 7.0 7.55 0:043 0.056 
sai 4.8 7.60 0.060 0.078 
#119 4.2 7.55 0.070 0.078 
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TABEE 1 (Cont "d) 


7 
Run ; : vA ; V 
(CPS) (psf) (VRMS) 
x 102 
jaa 3.6 5.98 0.055 0.0594 
ne Sal 5.98 0.043 0.0458 
4-3 7.9 5.92 0.030 0.0344 
4-4 10.1 5.92 0.040 0.0504 
4-5 12.6 5.92 0.030 0.0412 
4-6 14.0 5.92 0.019 0.0322 
a7 1726 5.88 0.016 0.0300 
4-8 19.3 5.98 0.016 0.0346 
4-9 15.6 5.92 0.017 0.0298 
4-10 lel 5.92 0.021 0.0344 
Pe 9.1 5.92 0.034 0.0392 
ane 6.6 5.92 0.033 0.0390 
ane 4.6 5.98 0.046 0.0504 
4-14 4.0 5.98 0.050 0.0574 
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TABLE 1 (Cont'd) 


x 
re f q en 2 : 
(CPS) (psf) (VRMS) 
x 102 
10-1 3.6 0.99 0.0039 0.0288 
to 2 5.5 0.99 0.0037 0.0266 
10-3 Ws 0.99 0.0038 0.0266 
10-4 9.7 0.99 0.0041 0.0308 
10-5 12.3 0.99 0.0039 0.0236 
10-6 14.5 0.99 0.0030 0.0206 
197 16.7 0.99 0.0035 0.0194 
10-8 19.2 0.99 0.0032 0.0184 
10-9 7G 0.99 0.0032 0.0190 
10-10 ae 0.99 0.0036 0.0224 
ema 10.9 0.99 0.0039 0.0288 
10-12 8.4 0.99 0.0038 0.0266 
10-13 6.3 0.99 0.0037 0.0266 
10-14 4.3 0.99 0.0039 0.0288 
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Run* 


5 


i656 


ee 


Igli=3 


9 


10 


Pie 


ie 


List 3 


11-14 


*The run number is followed by a coordination number which 


(CPS) 


ml 2 
13k 
1.5 
Pou 
m.3 
Iljswars: 


14.2 


TAT 
es 


4.4 


TABLE 2 
REDUCED DATA: MEDIUM VANES 


(psf) 


.56 
»45 
-40 
~45 
40 
mS 
-40 
aon 
op 
aD 
40 
-40 
245 


7 L 


mA 


(VRMS) 
x 102 
048 
U8 
103:7 
» 036 
20312 
7083 
2033 
029 
- (OSs 
nuS3 
.040 
.030 
036 


043 


to identify data as recorded on magnetic tape. 


40 


<l<> 


282 


ieee 


286 


w2iZ8 


ee 


ne02 


. 10 


. 166 


«2alu2 


184 


ree 


eee 


s2iio 


2e 


is used 


Run 


12-1 


i 2=2 


12-3 


12-4 


cae 


12-6 


(CPS) 


1Ur? 
mol 
5". 5 
Kor5 
or 2 
Lyre 
14.1 
1 20 
9.3 


6.8 


TABLE 2 (Cont'd) 


Pac 
ome 
2.24 
Zeal 
2.03 
2 05 
oo 
2.08 
2.08 
Z see 


2S 
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“A 
(VRMS) 
x 102 
0.090 
0.063 
0.047 
0.060 
0.047 
0.053 
0.050 
0.048 
0.050 
0.048 
0.050 
Om0ay 
0.050 


0.077 


<l<¢7 


. 346 


- Joe 


250 


274 


Src 


moO 


250 


. 206 


3 y' 


NN 


. 246 


296 


. 268 


. 380 


TABEEe2 (cont a 


v 
Run . : ; v 
(CPS) (psf) (VRMS ) 
x 102 
13-1 3.6 2.96 0.120 0.388 
lie? 5.6 2.65 0.077 0.376 
13-3 8.8 2.44 0.053 0.276 
13-4 12.0 2.65 0.077 0.332 
13-5 15.4 2.44 0.053 0.240 
13-6 18.2 2.39 0.060 0.290 
13-7 19.1 2.86 0.055 0.230 
13-8 16.9 2.34 0.060 0.282 
13-9 13.8 2.49 0.057 0.266 
13-10 10.8 2.65 0.080 0.342 
13-11 9.4 2.49 0.063 0.294 
13-12 Pel 2.44 0.053 0.262 
13-13 6.6 2.54 0.060 0.290 
13-14 4.5 70 0.100 0.375 
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TABLE 2 (Cont'd) 


: 
Run i : “A *y 
(CPS) (psf) (VRMS) 
x 102 
14-1 3.6 on73 0.015 0.204 
14-2 5.8 0.73 0.017 ' 0.204 
14-3 8.3 0.73 0.016 0.181 
14-4 118 0.68 0.022 0.219 
14-5 és, 1 0.68 0.016 0.142 
14-6 16.6 0.68 0.017 0.156 
14-7 20.5 0.68 0.015 OH 194 
14-8 18.3 0.68 0.016 0.148 
14-9 iSvs2 0.73 0.017 0.142 
14-10 2 7 lous 0.018 0.158 
14-11 9.7 0.03 0.017 0.188 
14-12 72 Og 0.016 0.188 
14-13 6.2 0.73 0.017 0.204 
14-14 5.3 0.73 0.016 0.204 
14-15 40 0.78 0.016 0.220 
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TABLE 2 (Cont'd) 


v 
Run ; p A *y 
(CPS) (psf) (VRMS ) 
x 102 
15=1 5.6 oe 20 0.140 0.428 
15-2 6.2 Ze 1. 0.090 ' " 0.398 
15-3 ono 270 0.598 0.246 
15-4 12.5 2.86 0.080 0.316 
15-5 Lom 2.76 0.060 0.256 
15-6 Py 2 2.54 0.070 0.294 
Loe 7 20.8 3me8 0.058 On282 
[3-8 ro. 5 2.60 0.070 Om292 
15-9 16.2 2.60 0.066 0.264 
15-10 io ae? 0 0.063 Og2i7'2 
15-11 ibs 3 2296 0.090 0i5383:6 
15-12 8.8 ZnO) 0.053 0.258 
15-13 Wied ZO 0.066 0.298 
15-14 4.9 Bol 0.120 0.390 
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Run* 


325 


5-6 


3-9 


5-8 


5-10 


5-11 


5-12 


3343 


5-14 


*The run number is followed by a coordination number which is used 


(CPS) 


5. 6 


10.0 
2.3 
14.3 
17.@ 


O. 3 


13.0 
%. 3 
6x6 
4.7 


3.6 


TABLE 3 


(psf) 


1.40 
1.71 
1.56 
1.45 
1.04 
2.08 
1.04 
1.56 
1.45 


1.66 


REDUCED DATA: LARGE VANES 


(VRMS) 
x 102 


0.150 


0.120 


0.080 


0.110 


0.100 


0.100 


0.120 


Or. 098 


0.120 


0.090 


0.079 


0.170 


0.130 


0.150 


to identify data as recorded on magnetic tape. 
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No 
<|1<7 


1.084 
0.976 
0 2922 
0.894 
0.876 
0.776 
t.OR2 
0.704 
0.960 
0.876 
0.776 
0.910 
lLuselultS 


1.108 


TABLE 3 (Cont'd) 


v 
Run ; : “A *y 
(CPS) (psf) (VRMS ) 
x 102 
6-1 3.6 mor 0.145 1.040 
6-2 Si T 36 ®:.'T20 1.186 
6-3 1 1.40 0.077 0.870 
6-4 10.0 Lal 0.100 0.842 
Ge5 12.4 1.48 0.101 0.908 
6-6 14.6 1.30 0.110 0.906 
G7 H6,.'8 1.04 Carez0 1... PSs 
6-8 Pot i, 0.090 0.584 
6-9 i> .2 1 a0 0.108 1.004 
6-20 i .6 1.40 0.090 0.816 
oe ml. 3 1 UO 0.928 
6-12 5. TZ 0.082 0. 9e2 
6=13 6.5 1.36 0.090 1 #082 
6-14 4.4 1.66 Oi7e20 L. 200 
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Run 


Ta 


7? 


1=3 


7-4 


Eee 


7-6 


7-7 


(CPS) 


10.3 
v3...6 
15.6 
17.7 
21.0 
16.2 
14.2 


eS 


(psf) 


1.25 
1.30 
1.14 
1.20 
1.04 
1.14 
1.04 
F.30 
1.45 


1.14 


TABLE 3 (Cont'd) 
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CA 


(VRMS ) 
x 102 


0.130 


0.100 


0.077 


0.083 


0.097 


0.097 


0.110 


0.083 


0.110 


0.083 


0.100 


0.081 


0.077 


0.110 


NO 
<|l<¢> 


0.924 


1.076 


0.850 


0.784 


0.766 


0.766 


0.986 


0.626 


0.876 


0.672 


0.832 


0.784 


0.816 


1 «062 


TABLE 3 (Cont'd) 


v 
Run : ; “A *y 
(CPS) (psf) (VRMS ) 
x 102 
8-1 4.0 1.61 0.109 0.864 
8-2 oy) 1.40 0.087 0 922 
8~3 8.4 re 0.066 0.708 
8-4 10.6 V.a4 0.083 0.724 
oO Eid ‘Laren 0.083 0.7652 
8-6 16.0 1 oer 0.083 0.638 
o=7/ mr. lew 0.077 0.480 
8-8 18.3 1.46 0.110 0.874 
8-9 14.1 er 0.077 0.596 
8-10 P23 T. 50 0.090 0.720 
8-11 ee Leer 0.073 0.682 
8-12 1 Teed 0.070 0.734 
8-13 Oye 1.30 0.090 0.920 
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TABLE 3 (Cont'd) 


v 
Run ; : “A ; v 
(CPS) (psf) (VRMS) 
x 102 
9-1 3.6 0.68 0.047 0.670 
9-2 6a 0.62 0.050 0.618 
9-3 9.6 0.57 0.063 0.702 
a 11.8 0.62 0.043 G.. “ae 
9-5 oy, 0.52 0.066 0.622 
9-6 18.1 0.47 0.073 0.674 
ony 20.9 0.52 0.044 0.424 
9-8 lem Geae 0.059 0.626 
9-9 13.4 0.62 0.039 0.388 
9-10 10.6 0.62 0.060 0.644 
oats 8.5 0.57 0.043 0.540 
9-12 eS 0.62 0.047 0.594 
daaa 5.5 0.68 0.050 0.670 
9-14 iO 0.68 0.043 0.670 
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q = 4.19 psf 
£ 
Run 
(CPS) 
1=1 3.6 
l=2 5.0 
13 la 
1-4 10 .@ 
=) 15270 
1-6 16.3 
i lore 
1=8 Lye 
= ioe 
1-10 i. 6 
dye]. Dien? 
Tar2 6.4 
1-13 4.2 
Peis 3150 


TABLE 4 
Experimental Phase Difference Between 
Anemometer and Hot Wire 


Small Vanes 


DB if 


Run 

(degrees) (CPS) 
7 2-1 37.0 
lz 2-2 6.0 
ils 2-3 8.8 
20 2-4 ale, 2 
iy 2=5 io. 3 
22 2-6 15.0 
27 2-7 16.9 
22 2-8 eo, 1 
23 2-9 18.1 
18 2-10 a9 
16 2-11 14.1 
iP 2-12 on 
9 2-13 7-0 
6 2-14 5.0 


50 


q = 2.34 psf 


P 


(degrees) 


aa) 
13 
ak 
2 
Di2 
24 
24 
Zi] 
29 
28 
22 
12 
Is 


12 


TABLE 4 (Cont'd) 


de=) 7.509 pet q = 50 psf 
Run . P Run : @ 
(CPS) (degrees) (CPS) (degrees) 
3-1 S60 11 4-1 340 5 
3-2 Jac i 4-2 Sia 7 
3-3 8.6 Zul 4-3 Jig 8 
3-4 Pree ao 4-4 LO 10 
3-25 L3G 32 45 2 .@ 18 
3-6 1 36 4-6 14.0 7 
3-7 LOmZ 312 4-7 le? 0 Zl 
3-8 a. 7 ay 4-8 19.3 iL 
3-9 12.4 oy 4-9 L5 Zi 
3-10 Dry 23 4-10 1a 19 
3-11 70 20 4-11 Sal. 9 
3-12 4.8 15 4-12 6.6 8 
3-13 4.2 ins 4-13 4.6 ms 
4-14 4.0 4, 
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Run 


10-2 


10-2 


10-3 


10-4 


i= 


10-6 


10-7 


10-8 


10-9 


10-10 


10-11 


10-2 


10-13 


10-14 


TABLE 4 (Cont'd) 


q = 1.00 psf 
g 
(CPS) (degrees) 
326 23 
Sie 29 
7.4 so 
og a 
hog 36 
14.5 38 
F6e7 49 
or 56 
17.6 48 
1375 ov 
10.3 34 
8.4 30 
65 29 
Lt 26 
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qy= 1.43 psr 
f 
Run 
(GPRS) 
= 1 2.7 
L1l-2 6.4 
ea 8.6 
11-4 12 
li-5 esl 
11-6 5.5 
ie ee). 2 
lt-3 193 
11-9 16.8 
11-10 14.2 
11-11 G7 
11-12 a7 
11-13 Seo 
11-14 4.4 


Experimental Phase Difference Between 


Anemometer and Hot Wire 


D 


TABLE 


5 


Medium Vanes 


(degrees) 


ik 


46 


30 


27 


30 


oi 


Run 


12-1 
12e2 
12-3 


12-4 


12-10 


53 


12-11 


12-12 


Tae 


12-14 


CGPS) 


ee | 
15.5 
re.5 
9), 2 
Vee 
14.1 


i a8 


q. S222 3epst 


Pp 


(degrees) 


40 
23 
24 
22 


ze 


q "=" 2, 50mpst 
f£ 
Run 
(CPS) 
oe 2.6 
eZ Sal 
133 8.8 
13-4 Pe, 0 
1325 15.4 
13=6 Ve. 2 
ie hy IE 
Ea re. 
ife9 13.8 
1310 EUS 
13-11 9.4 
13-12 weet 
13-13 oO. 0 
13-14 4.5 


D 


(degrees) 


18 


28 


22 


36 


TABLE 5 (Cont'd) 


Run 


14-4 
1445 
ee 


14-7 


Lael 
14-12 
14-13 
14-14 


14-15 


a. 3 


14.1 


16.6 


20,5 


re. 3 


Ee, 2 


re. / 


q = 0.71 psf 


P 


(degrees) 


16 


a5 


20 
42 
34 
igh 


44 


21 
12 
13 
5 


14 


TABLE 5 (Cont'd) 








q = 2.85 psf 
Run ; P 

(CPS) (degrees) 
15-1 2 E 
15-2 OZ 2 
15-3 9.9 0 
15-4 B25 0 
15-5 i. 1 4 
15-6 eye 40 
15-7 2066 20 
15-8 To.5 27 
15-9 Le 2 33 
15-10 5.8 0 
Po Pt Laie 0 
15-12 8.8 1 
15-13 icles. 0 
15-14 4.9 0 


TABLE 6 


COMPARISON OF FOURIER COMPONENTS 


Run S ; Signal Fundamental 2nd Harmonic 3rd Harmonic 
(x 10°) Dy Ayla, Qo = Ag/A, Qs 
(degrees) (%) (degrees) (%) (degrees) 


1-1 Stee? ey = m7 - 4.4 - 
en -2.0 8.5 -13 TORY 47 
1-2 4.83 ey = 28.7 - 3.4 - 
eA V8 30.8 56 2.8 ZS 
1-2* 4.83 ey = 50.5 : 6.0 - 
Qn Vad 24.6 20 9.4 -69 
t=3 7.45 ey - 27.4 - Ze - 
en 14.0 Or. 5 6.4 138 
1-4 962 ey - 10.0 - 2.8 - 
en 7.0 2.0 -104 1.6 162 
Note: 
i P, is phase difference between hot wire and anemometer 


Signal for component in subscripts (positive for hot 
wire leading anemometer). 


Z a ) indicates amplitude ratio of components in subscripts 
” ) for specified signal, i.e., e, or Qa: 


* indicates different waveform sample in run. 
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TABLE 6 (Cont'd) 


Run S Signal Fundamental 2nd Harmonic 3rd Harmonic 
(x 102) A A A,/A, Dz 
(degrees) (4) (degrees) (%) (degrees) 
2-1 4.69 ey - 2009 - ei - 
on 4 Feso IO Sade 56 
2-14 Oa ey - Diu 2 = 4. - 
en 14 mi.oe -13 9 4 Sab 
2-13 6.51 en - Zoe - 2.8 7 
en 28 8.4 50 ae 5 
Za 8 lZeZ5 ey - ZC a - 4.4 S 
CA OFZ 8.8 D7 7.8 -16 
2-4 1 SS en = Tony - Zao - 
en 7 ho 2 lesa feee 
2-4 14.55 ey - Ge - Bac = 
CA 4 toes) 25 | ee 
2-11 Fore en ~ 24.5 = ag - 
on 12 Jo 1 205 7 


ey! 


Run 


16-8 


16-9 


16-10 


16-11 


16-12 


16=13 


16-14 


He=t5 


16-16 


Lé=17 


16-18 


16-19 


16-20 


TABLE 7 


REDUCED DATA: SMALL VANES 


£ q en 
(CPS) (psf) (VRMS) 
x 10 
36 dus 2 0.0310 
5.4 4.26 0.0240 
“year 4.26 0.0200 
8.4 4.26 0.0170 
9.4 Ly te 0.0180 
10.8 4.26 0.0200 
bis 4.26 0.0140 
£36 4.26 0.0104 
14.2 4.26 0.0104 
15.6 4.26 0.0104 
16.6 4.26 0.0096 
7 7 4.26 0.0097 
19%) ZO 0.0095 
18.3 4.26 0.0096 
16.3 4.26 0.0099 
137-7 4.26 0.0104 
Tee 4.26 0.0148 
29 4.26 0.0191 
7.4 4.26 0.0180 
rl 4.26 Onoz12 
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<l<> 


0/728 


0562 


0474 


0444 


0526 


0621 


. 08Us 


0416 


Oey 


0887 


20507 


US72 


0444 


0411 


0387 


0406 


0511 


. UD 


0454 


0507 


Run 


De /igeel 
ies 
1jes 


17-4 


Pew 
17-11 
1?=12 
ce13 
Lyme 
bred 
l= 16 
ey 


Lie 8 


(CPS) 


3.6 


6.8 
onl 
95 


10.4 


$3.1 
14.2 
5.7 
3.5 
pO 1 
Is .0 


6.6 


TABLE 7 (Cont'd) 


(psf) 


6.03 
508 
5.28 


5298 


5.98 
5 0S 
5 08 
rs 
BZ 
D 76 
5 6 
5.6 
By i 2 
5S 
Divo 
».O8 


2» 106 


59 


fA 


(VRMS) 
x 102 


0.0524 
0) e592 
0.0320 
O 0279 
0.0341 
0.0382 
0.0247 
0. 0183 
0.0150 
030139 
O2Ciy 
0.0139 
0 S019 
0.0139 
0.0183 
0.0341 
0 0289 


0.0310 


INO 
<i<> 


OisO793 
0.0609 
0.0507 
0.0473 
0.0588 
0.0684 
0. 0527 
0.0456 
0.0429 
0.0411 
0.0470 
0.0489 
0.0454 
0.0428 
0.0451] 
0.0647 
0.0533 


0.0505 


Run 


18-1 


18-2 


ie 3 


n= 10 


eee) 


18-12 


hee 13 


18-14 


bed 5 


18-16 


18-17 


18-18 


Lae l9 


18-20 


(CPS) 


Sb 


10.7 
11.4 
eye 2 
12.8 
Lome 
14.0 
5.7 
7.2 
bOml 
18.3 
15.5 


Lien/ 


845 


(psf) 


73 
Peo 1 
a1 
Ta 3 
Tap 3 
7.48 
7.48 
7.48 
7.48 
73 
7.48 
7.48 
7.00 
lao 
7.48 


7.48 


TABLE 7 (Cont'd) 
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oh 


(VRMS ) 
x 102 


0.0784 
0.0469 
0.0428 
0.0398 
0.0377 
0.0449 
0.0531 
0.0459 
0.0346 
0.0314 
0". 0283 
0.0251 
020173 
0.0162 
0.07 3 
M062 
Unut7 3 
0.0251 
0.0490 


020567 


NO 
<ij<d> 


0.0836 
0.0562 
0.0540 
0.0499 
0.0507 
0.0584 
0.0709 
0.0640 
0.0546 
0.0499 
0.0491 
0.0454 
0.0449 
0.0465 
0.0554 
0.0502 
0.0445 
0.0463 
0.0669 


0.0497 





Run 


19-4 


19=5 


19-6 


19-7 


ig=s 


19-9 


i= 


19-11 


Een 2 


iSe15 


19-14 


E915 


19-16 


19-17 


Tos 


Pao 


(CPS) 


36 
5.5 
7a 
8.6 
9.6 
Ledeen 
bee 7 


14.2 


18.0 
17.2 
15 
13.8 


12.0 


TABLE 8 


2.08 


2.08 


eS 


2.08 


2. 


2.08 


mS 


2.08 


2.06 


Zao 
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REDUCED DATA: MEDIUM VANES 


SA 


(VRMS 
sli) 


0.084 
0.058 
0.043 
0.041 
0.052 
GEU55 
0.044 
0.046 
0.050 
0.046 
0.040 
0.046 
0.047 
0.049 
0.044 
0.043 
0.048 
0.060 


0.040 


N 
<l<> 


©. 370 


0.326 


OPZ36 


0. 280 


0.266 


0.280 


0.222 


0.216 


0.246 


0.236 


0.206 


@.232 


0. 288 


0.236 


0.214 


0.216 


0.248 


Omol2 


W226 


Run 


20-10 


20-11 


20-12 


20-13 


20-14 


20-15 


(CPS) 


6a 


ga 


10.6 


15.0 
16.4 


18.4 


14.1 
(ea 


8.0 


(psf) 


2. 70 
281 
3402 
2.86 
203 
25/0 
2.70 
eel 7 


3.54 


200 


2./0 


TABLE 8 (Cont'd) 
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CA 


(VRMS) 
x 102 


0.138 


0.064 


0.056 


e057 


0.091 


0.062 


0.057 


0.065 


0.066 


0.061 


0.057 


0.064 


02059 


02072 


0.054 


ie) 
<ij<¢> 


0.440 
0.270 
0.248 
0. 240 
0.352 
0.258 
0.236 
0.254 
0.276 
0.258 
0.238 
0.270 
0.234 
0.300 


0.238 
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FIG. 2: ANEMOMETER STATIC CALIBRATION 
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FIG. 3: ANEMOMETER INERTIAL RESPONSE 
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FIG. 4: SHUTTER VANE BLOCKAGE 
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FIG.5: DIAGRAM OF DATA ACQUISITION SYSTEM 
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PIG 7: SENSOR INSTACEATION. DETAILS 
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PRESSURE (W.C.R. TUNNEL) 
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FIG.9: MEAN ANEMOMETER DRAG vs. MEAN 


DYNAMIC PRESSURE (NONSTEADY TUNNEL) 
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FIG.10: UNSTEADY “DRAG VS. FREQUENCY 
(SMALL VANES) 
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FIG.12: UNSTEADY DRAG vs. FREQUENCY 
(LARGE VANES) 
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FIG.13: ANEMOMETER PHASE LAG (SMALL VANES) 
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Sweep Rate = 20 msec/cm 





Sweep Rate 29 msec/cm 
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FIG.18: DRAG RATIO VARIATION (Small Vanes, q=1.0 psf) 
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FIG.19: DRAG RATIO VARIATION (Small Vanes, q=23 psf ) 
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FIG. 20: DRAG RATIO VARIATION (Small Vanes, q=4.2 psf) 
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FIG.21: DRAG RATIO VARIATION (Small Vanes,q=5.9 psf) 


83 






0.10 bs UNTAGGED SYM-RUN 3 


TAGGED " =- 18 
© ~ ANEMOMETER 


008 xX Oo A~QUASI-STEADY 
BR. DRE 

0.06 

0.04 

0.02 

NOTE SUPPRESSED 
O 
0.02 
4 8 “4 1 20 


Frequency (cps) 





FIG.22: DRAG RATIO VARIATION (Small Vanes, q=7.6psf) 
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FIG.23: DRAG RATIO VARIATION (Medium Vanes,q=0.7 psf) 
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FIG. 24: DRAG RATIO VARIATION (Medium Vanes,q=1,4 psf) 
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FIG.25: DRAG RATIO VARIATION (Medium Vanes, q=2.2 psf) 
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FIG.26: DRAG RATIO VARIATION (Medium Vanes, q=2 6pstf) 
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FIG.27: DRAG RATIO VARIATION (Medium Vanes, q=2 .8psf) 
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FIG.28: DRAG RATIO VARIATION (Large Vanes,q=0.6 psf) 
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FIG. 29: DRAG RATIO VARIATION (Large Vanes,q=1.2 psf) 
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FIG.30: DRAG RATIO VARIATION (Large Vanes, q=1.3 psf) 
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FIG.31: DRAG RATIO VARIATION (Large Vanes,q=1.5psf) 
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FIG.32: DRAG RATIO VARIATION (Large Vanes, q=1.7psf ) 
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FIG.35: STATHAM ACCELEROMETER CALIBRATION 
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APPENDIX A 
Accelerometer Calibration: 

Calibration of the Statham accelerometer was obtained by measur- 
ing amplitude and frequency of motion, plus phase shift, between 
motion and accelerometer output when the accelerometer was mounted 
on the table of a Calidyne shaker unit, model A88. 

Amplitude of motion was obtained using a Bently distance detector 
unit which is effectively a reluctance gage that senses the change in 
distance of the air gap between a flat, non-magnetic conducting mater- 
ial and the sensor coil face. For the calibration range segelitiead? 
the reluctance gage frequency response was the same as the static 
sensitivity with a distance resolution of + 20 micro-inches. Frequency 
was sensed to an accuracy of + 0.1 cycles per second (cps) using a 
conventional electronic counter unit, Hewlett~Packard, model 522. 
Phase angle difference was sensed by an AD-YU, type 405 phase meter 
which had an accuracy of + 1.0 degree. 

The Statham accelerometer was operated at a 10.0 volt (D.C.) 
bridge excitation level and had an approximate static sensitivity of 
4.0 millivolts per "g"' where the gravitational unit "g'' corresponds 
to a standard value of 32.2 feet per second2. The dynamic response 
of the accelerometer is shown on Figure 35 including both gain, Go Cw ); 
and phase lag characteristic Aga. The acceleration response was 
checked for linearity at constant frequency up to two g's. The un- 
damped resonant frequency, fy, was evaluated as being 218 cps using 
the criteria of ninety degrees phase lag. The accelerometer was 
estimated to have a damping ratio of 0.558 relative to critical damp- 


ing, and comparison is shown on the frequency response plot. 
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The calibration may be considered as related to true acceleration 
by the relationship: 
If true acceleration is 
X(t) = Asin wt 
Then the accelerometer reading will be (A.O1) 


e.(t) = AG. sin(wt-a¢,) 


Where G, and Aq g are obtained from Figure 32. 
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APPENDIX B 


Anemometer Corrections For Tunnel Floor Accelerations 


It is the purpose of this section to illustrate a method for cor- 
recting the drag sphere anemometer balance-voltage reading for the 
influence of tunnel floor accelerations in the horizontal axial di- 


rection. A sketch of the environment is indicated below: 


CLL LA pf LIDS IG TL LAL LIT eee 
















Reference | . 
Inlet Velocity | 
Station T 
Mounted 1 
V+ v(t) Anemometer 9 
A 
e — > ; 
Rotating 
Tunnel Shutter 
Supports Accelerometer Vewmet 


CON Ltt “UMbEteT 


Sketch of 2.0 x 240 foot test seetion 


Nonsteady Wind Tunnel (Not to Scale) 
At the reference velocity station, the hot wire provides an indi- 


cation of the air velocity. This station will be used as a source for 


the origin of time. 

V+v@)=V ri +O Sin wt | (B.01) 
where O© represents the magnitude of velocity perturbation relative to 
the steady velocity value. The floor will experience an acceleration 
due to the aerodynamic loading induced by the operation of the rotating 
shutter mechanism, and this will occur principally in the horizontal 
axial direction. For purposes of this discussion, let us assume that 


the motion is harmonic in character at the dominant forcing frequency. 
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The floor acceleration can be expressed as: 

x we) = C sin (wt - ¢,) (B.02) 
where it is recognized that the quantities Ci> and Pas are not 
measurable directly. 

The Statham accelerometer, which is firmly attached to the floor 
beneath the drag sphere anemometer installation, was aligned to re- 
spond to accelerations in the horizontal axial direction. Its elect- 
rical signal output due to the acceleration environment of equation 
B.02 will be: 

é.(t) = C, sin(wt-4) (B.03) 
where the quantities Cy and D> are directly measurable. From labora- 
tory shaker table calibrations, it is known that: 

ee = 9 G. (w) » and 

AQ. F h, - 

Also we know from a laboratory shaker table calibration of the 
drag sphere anemometer that 

Ae, (t) = C, sin(wt-4,) 

C; sin|(wt -&)- ag, | 
where = = C,G,(w) ) and Ad, = P,-® 


Combining equations B.02 through B.04 allows us to estimate an inertial 


(B.04) 


correction to the drag sphere balance reading. Note that Pi is can- 


celled out and is thus never required. 


Ae,(t) = C, aes sin] (wt 5 >) -A®, + a¢, | (B.05) 


The observed drag sphere anemometer balance output is directly 


measured and may be expressed as: 


= Bsin Con ~ $, ) (B.06) 


A oss 


oS 


and it is desired to obtain a corrected reading in the form: 


C nese _ nas Sin (wt _ y) 
= Cas ae 


A ops 


(B.07) 


Solving equation B.07 in terms of quantities involved in equations 


B.05 and B.06, one obtains: 


2 fe 
Bosee =| B'+C GP 28C, Zt cos(4,- 4-24, +04) 


and 


GY 
> 


sin($,+A,- Ag) | 


sin d, - 


a 


Y = arctan (8.68) 


G)} 6) 
> 





o}O) 0). 


cos  - 


cos( ¢,+ a4, -A¢, ) 


G) 


It should be noted that the phase relationship between the un- 
Steady velocity and drag (hot wire and anemometer) is the correct 
relationship between the velocity at the hot wire and the drag at the 
anemometer. However it is not the relationship between the velocity 
and drag at the same point. Due to the longitudinal separation of 
hot wire and anemometer the pressure signals producing velocity per- 
turbations reach the anemometer prior to reaching tie hot wire. Based 
on the speed of sound in air, the mean flow velocity in the tunnel, 
the frequency of oscillation, and the distance between hot wire and 
anemometer, the maximum phase difference between velocity signals at 
the hot wire and anemometer, for the range of test ee used, 
would be 7 degrees. This phase difference, due to the above factors, 
means that the phase angle (or time lag) between a velocity perturbaticr 


at the anemometer and the anemometer drag response is greater than 


that indicated by phase angle measurements between the hot wire and 
the anemometer. Figure 36 shows the computed phase shift, due to the 
axial separation of the hot wire and anemometer sensors,as a function 


of frequency for four representative values of dynamic pressure. 
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APPENDIX C 

Hotwire Response to Base Acceleration 

An elementary beam theory analysis has been made for the canti- 
levered rod (supporting the hot wire probe and containing electrical 
leads) to determine the velocity of the probe induced by inertial 
loads during tunnel floor vibrations in the horizontal plane. The 
entire length of the rod (shown in the sketch below) was assumed to 
be identical to the hollow stainless steel portion in specific weight, 


cross section, and elastic properties. 


& a Probe 
Epoxy Tip 
Steel Rid 


Tunnel Floor 











o@s2"G)) 


— | =e 
YW Olid: ey 


Rod Cross Section 


SKETCH OF HOTWIRE INSTALLATION 


The length of the rod was taken to be 6.5 inches since the rod 
is fixed at the bottom of the tunnel floor by a threaded collar 
arrangement. This is a conservative value because the floor pro- 
vides some support due to the close fit between the rod and the hole 


through the floor. 
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The following properties and values were chosen or computed 
for the analysis: 
Modulus of elasticity E = 29 x 10° lbs/in2 


0.284 lbs/in? 


Specific Weight x 


Second moment of the area with respect to the ZZ axis; 


J EL wy 1 
Iyz = © (dy*-dy") 
64 
Izz = 1.895 x 10 int 
Cross sectional area; 
Av= 7 (a, *#ané 
r (dy?~dy”) 
A = 0.0396 in? 
Length of rod; f = 6.5 in 
Weight of rod per unit length Ma =VA 


4A°= 0.01125 1bs/in 
Weight of electrical leads Lh = 0.05 10° 
Effective Weight MA = Wr’ + 4” 
AA = 0.0118 lbs/in 
The vertically mounted rod, subjected to an acceleration of one "g" 


in a horizontal direction, is equivalent to the rod mounted horizon- 


tally, loaded by its own weight, as shown in the following sketch. 


uniform load of 4£& Ilbs/in 
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Elementary beam theory gives the deflection of the free end 


(y = f ) as: 
3 
x = Wd (C.01) 





where, W = w4 under a one "os" condition, and the maximum deflection 


is then 


ws. = (040767) be. (6.5) ime 


° 8 (2.9 x 106) lbs/in2 (1.895 x 1074) in4 


0.00048 in 


=o 
If the load were assumed to be applied harmonically, ~i.e., 

W= Wo, sinwt, (C.02) 
then the deflections of the free end can be approximately expressed 
by 

xo = 2oeesan ae (C.03) 
provided the frequencies of interest are small compared to the lowest 


natural frequency of the rod. 


The natural frequencies of a continuous uniform beam are given by 


2 
aa (G4) | faa (C..04) 


where for a cantilever, BA = 1.875 is the lowest eigenvalue. 





Mass per unit length 
m = w/g 


3.06 x 107? lb sec?/in2 


m 


Substitution of values into equation C.04 yields 


29 x 10© 1bs/in? (1.895 x 10%) in‘ 
(1.875)? (3.06 x 1072)1b sec4/in4 (6.5)4 in4 


CG), = 
J, = 1113 radians per sec 
or £, = 177 cycles per second 
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Since the maximum frequency of interest is 20 cycles per second, 
the ratio (f/f) is sufficiently small for an order of magnitude 
analysis based on (C.03): 


For x = X) sin wt, the velocity of the probe tip is x =Wx, cos at. 


O 

Thus the tip velocity, relative to the base, can be further ex- 
pressed as x = ndxX, cos wt where n is the number of "o's" acceler- 
ation. 

Converting units from radians to cycles and inches to feet yields 
an expression for the maximum tip velocity in terms of the test varia- 
bles n ("g" factor) and f (frequency): Bi ase = 6.36 nf x 107° ft/sec 

During the nonsteady tunnel test program the maximum horizontal 


acceleration recorded was one-fourth "g" (n = 1/4), while maximum 


frequencies were 20 cycles per second; thus 


Kmax = 6+36 (1/4) (20) x 107° £t/sec 


-6 
a. 31.8 x 10 ft/sec 


The minimum velocity perturbations which occur anywhere in the 


investigation are about 0.26 feet per second. Therefore : = repre- 


- 


sents about 0.01 percent and is negligible 
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APPENDIX D 


Tunnel Floor Velocity Eifects on Hotwire Signal 


The following analysis was used to determine the velocity pro- 
duced at the hot wire as a result of tunnel floor velocity, considering 
the hot wire support to be rigid. 

The floor motion was assumed to be simple harmonic, and expressed 
as 

~ = Asin wt (D.01) 
from which the floor velocity and acceleration are given by 


ee 
i! 


and 


: 2 : 
X=-w Asin wt (D.03) 
respectively. 


24 can be determined from the accelerometer data 


The factor WJ 
for each run. Then using the oscillation frequency of that run the 
maximum displacement, A, can be determined, and the product A is 
the velocity amplitude. Due to the rather large variations in accel- 
eration, and the fact that the term A varies with frequency, no ex- 
pression can be obtained which describes the entire test. Instead 
the quantity A must be computed for each test point. 


The following sample calculation indicates the method used: 


Run 5-1: 


from raw data: f 3.7 cycles per second 


0.2 millivolts (RMS) 


ec 


gaecelerom@ter sensitivity: k 9~>“ajazemillivolts per “g” 


Coe hi= Ce). = 0.0498 "g's"! (RMS) 


2 
wA = 1.60 feet/sec2 (RMS) 
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Now the velocity amplitude, wA, is given by 


XY =wAs WA _ 160 
F Gd) ent 


X cemaldiOS8 Sot secu hi®MS) 
From previous data reduction 


A 
€,,= 0.21 volts (RMS) => v=24 fthec (RMS) 
The relative magnitude is then 
Xe /s = O.@0B6 7. 
or about 0.3 percent. 
The maximum value of %p/V found in the test runs was of the 


order of 2 percent, which was considered to be a negligible influence. 
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APPENDIX E 
Quasi-Steady Theory: 

The ability of the drag sphere to sense dynamic pressure changes 
by the change in drag force is clear for "steady" velocity changes, 
however, the ability of the time-varying drag force to provide a 
reasonable representation of time-dependent variation of dynamic 
pressure cannot be assured by a linsambas approach unless verified 
by experiment. The quasi-steady theory, as explained here, is based 
upon the premise that the drag force coefficient remains invariant 
and that the instantaneous value of drag is dependent directly upon 
the instantaneous value of dynamic pressure. Comparison between exper- 
iment and quasi-steady theory is the purpose of this research, with 
the goal of showing both the drag magnitude and phase shift as a 
function of dimensionless frequency. Quasi-steady theory assumes no 
phase shift. 

If the velocity is assumed to be a harmonic component superim- 
posed upon a steady term, then we have that: 

V+v(t) = V(1+ osinwt) (E.01) 
where @ represents the amplitude ratio between unsteady and steady 
velocity. The instantaneous value of dynamic pressure can then be 


expressed as: 


% sev “| (14 = =) +2o0 sin wt - x COS 2ut| (E.02) 


Applying the quasi-steady assumption to drag force yields: 


D+ d(t) aa AV ‘|(+ F)+20 Sin wt - F cos 2ut | (E.03) 
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It will be noted that the steady drag term increases slightly over 
the value in uniform flow due to the presence of the 6472 term in 
equation E.03 . As explained in Section 3, the influence of the term 
was experimentally verified as being small, which is reasonable to 
expect for velocity perturbations up to about 15 percent of the steady 
veloeity peiwe, , | v(t S O15 

From equation E.03, we can readily discern that the ratio of 


unsteady drag to steady drag by quasi-steady theory is given by: 


2 
a =2osin wt - = cos Cut cemran) 


Next we consider equation E,04 from the standpoint of mean square 
values, which becomes upon applying the limit operation for the mean 


Square of a temporally varying term: 


‘ 
= bimet +/ (A c2sintet -2osinet cos 2ut + 





re 
D = T+ 
T cos* 2et) dt } 
dé _ 2 ios 
= Pot(i+ & (E.05) 


U 


a” 


Defining d =Jd° and considering equation E.05, 


4 V, 
LeFo(i+F =VZo(i+ 5+) 


We now must recognize that the velocity perturbation © corresponds 
to a peak amplitude and may be expressed in terms of the RMS of 
velocity perturbataon, — , by OF =/2 € : 
Hence we obtain that 

£2e(i+ Se.) 
which to a first order term may be stated as; 
=2 


Equation E.06 is a statement of the quasi-steady drag ratio as used 


(E.06) 


tylaQ> 
<|< > 


in Section 3 of this analysis. 
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APPENDIX F 


Fourier Analysis of Signal: 


A simple analysis was made of the time varying electrical signals 
from the hot wire and drag sphere anemometer in order to evaluate the 
quality of the unsteady velocity perturbation and the response of the 
anemometer to the various harmonics. Strictly speaking, the analysis 
using a Fourier series of consecutive harmonics is not valid for this 
experimental program, since it is possible that a frequency may be 
present in the signal which is not an integer multiple of the funda- 
mental, and hence a Fourier analysis using a truncated series could 
give a false indication of signal content. However, the Fourier 
approach was rationalized in this instance on the basis that the time 
to do a true power spectral and cross-spectral analysis was not con- 
sistent with the available time, and at least an approximate indication 
of frequency content had to be obtained. 

Consider a periodic function of period 2T as shown in the sketch 


below: 


+(x) 


Zul, 
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The Fourier series representation for the cyclical function may 


be expressed as: 


£(x) 











oo 
Pe + > (A. cos OT +4 By sin nrX (F.01) 


2. ee i Ts 
where , et 
ee x ™ 
Ay = + { f(x) )cos — lx poy Sali, Beets (F202) 
oer 
B, = — (¢ $1) sin SES dx pee ven: Uanas (F.03) 
One complete cycle of the saat to be analyzed was obtained 


from an oscilloscope photograph (cf., Figure 30) and divided into 


36 equal parts, such that 


ul 


36 AX; CT radians = 360 degrees 


LK, = 10 degrees 
X; = 10i degrees (i= 0,1,2,-"',35) 


t 


Now express period as 


2T=360 degrees 
Then at each 10 degree station (x; ) through the cycle, the amplitude, 
£(x;); wae measured. This information allowed the coefficients, 
A,, and B,, to be evaluated, using a numerical approximation to the 


integral te * - and F.03, as 


== es F(X;) cos NX; d%; (F.04) 
(20 
which, with x; “w degrees, is 
A. ie £(x:) cos NX; (F.05) 
l= 
Similarly 
LS f(a) si 
-— Xi) sin nx; (F.06) 
By 18 io 
Equation (F.01) may be alternately expressed as 
_A..S 
f(x)= Se + 2 C, cos ( w,t - 4,) (F.07) 
where 
z 2 oe 
C. =/A, +B) (F.08) 
and B 
¢ = arctan EE (F.09) 
n 


ISLS 


The constant term, Aj, is a function of the axis location (from 
which £(x,) is measured) and is of no value in this analysis. Only 
three components were computed, i.e., the fundamental, second, and 
third harmonic. After C, was computed (n = 1,2,3) for a given signal, 
the magnitudes of the second and third harmonic components, relative 
to the fundamental, was given by the ratios, Cy/C, and C3/C,. Compu- 
tation of dq, provided phase relationships among components of one 
signal, and also phase relationships between two signals processed 


relative to the same time scale and origin. 
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